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Coulomb coupling between proximal layers in graphene heterostructures results in efficient energy 
transfer between the layers. We predict that, in the presence of correlated density inhomogeneities 
in the layers, vertical energy transfer has a strong impact on lateral charge transport. In particular, 
for Coulomb drag it dominates over the conventional momentum drag near zero doping. The 
dependence on doping and temperature, which is different for the two drag mechanisms, can be 
used to separate these mechanisms in experiment. We predict distinct features such as a peak at 
zero doping and a multiple sign reversal, which provide diagnostics for this new drag mechanism. 



Recently developed vertical heterostructures [T] com- 
prised of a few graphene layers separated by an atom- 
ically thin insulating layer afford new ways to probe the 
effects of electron interactions at the nanoscale. Typ- 
ical layer separation in these structures (1-2 nm) can 
be very small compared to the characteristic electron 
lengthscales such as the de Broglie wavelength and the 
screening length. This defines a new strong-coupling 
regime wherein the interlayer and intralayer interactions 
are almost equally strong. Fast momentum transfer be- 
tween electron subsystems in the two layers and strong 
Coulomb drag have been predicted in this regime [2—7; 
with characteristic dependence on doping, temperature 
and layer separation distinct from that in previously 
studied systems [8]. Recently, measurements of strong 
drag in graphene were reported |H [TO]- 

In this article we focus on another effect that becomes 
prominent in the strong coupling regime: vertical energy 
transfer mediated by interlayer electron-electron scatter- 
ing. We predict that this process can give rise to lat- 
eral energy flow in the electron system, which under the 
conditions discussed below can directly impact electric 
transport. In particular, it leads to a characteristic con- 
tribution to the in-plane resistivity and dominates drag 
near charge neutrality (CN), see Fig{l]and Figj3j 

Interlayer electron-electron scattering, which governs 
drag, is strongest at CN since the long-range Coulomb 
interactions become unscreened near CN. The scattering 
rate decreases with doping away from CN as 7 oc t/ _1 (/i), 
where is the massless Dirac density of states, /1 3> 
ksT is the chemical potential, see Eq.([8|. The 1/v scal- 
ing is completely analogous to that found for intralayer 
scattering jTTJ [T2] . Crucially, while both the energy and 
momentum transfer rates peak at CN, their impact on 
drag is markedly different. Since the sign of momentum 
drag depends on the polarity of charge carriers mo- 
mentum drag vanishes at CN [2HZ]. In contrast, energy- 
driven drag features a peak at CN, see Fig(3] 

Our drag mechanism arises due to coupling between 
vertical energy transfer and lateral charge and energy 
transport via spatial density inhomogeneity which is in- 
trinsic to graphene. Density inhomogeneity is known to 




FIG. 1: Different mechanisms for Coulomb drag in graphene 
heterostructures. Energy-driven drag dominates over momen- 
tum drag near zero doping, whereas momentum drag domi- 
nates at higher doping. The sign of the drag response depends 
on carrier polarity (a). For potential fluctuations of equal sign 
in the two layers, Eq.|l]), the net drag (b) features a pair of 
nodal lines with an avoided crossing. The drag sign at zero 
doping is positive, as in the bipolar regime. This is distinct 
from momentum-only contribution to drag (c) smeared by 
correlated density fuctuations, <5/ii w 8^2 . 

be particularly strong at CN in the electron-hole pud- 
dle regime [H]. providing the dominant disorder poten- 
tial in clean samples. When a charge current is applied in 
layer 1, density inhomogeneity produces spatially varying 
heating/cooling [see Eq.|2])]. Strong thermal coupling be- 
tween the electron systems in the two layers, mediated by 
the interlayer energy transfer, leads to a temperature pat- 
tern in layer 2 that tracks that in layer 1, ST2(r) ps 5T\(r). 
Further, since the disorder correlation length £dis can 
reach 100 nm in G/BN heterostructures [T5l fT6] . exceed- 
ing the layer separation by orders of magnitude, the po- 
tential fluctuations are nearly identical in the two layers, 

(6^(r)SfM 2 (r')) > (1) 

for r w r' . As a result, the position-dependent ther- 
mopower induced by the gradient V5T2{r) is correlated 
with the heating/cooling pattern in layer 1, giving rise to 
a nonzero ensemble- averaged drag voltage in layer 2. 

Our mechanism predicts a particular sign of the energy 
contribution to drag. As a result, the density dependence 
for the net drag (energy and momentum combined) fea- 
tures a split-up pattern of nodal lines with an "avoided 
crossing" at zero doping, as illustrated in Fig. [I] (b) . The 
double sign change along the main diagonal n\ = 712 and 
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FIG. 2: Feynman diagrams for momentum (a) and energy 
(b) contributions to drag. Wavy lines represent interactions, 
dashed line represents disorder averaging. The ladder in (b) 
represents a long-wavelength charge-neutral mode. 

the peak at rix,2 = make the energy-driven drag easy 
to distinguish experimentally. 

As a parenthetical remark, the correlated density in- 
homogeneity, Eq. ([!]), also affects the momentum drag, 
however its effect is opposite to that of the energy con- 
tribution (see FigJI] b and c). If momentum drag were 
the dominant contribution near zero doping, the pattern 
of nodal lines would be such that the drag sign was con- 
stant along the main diagonal. This qualitative difference 
makes it easy to differentiate between the two cases. 

Below we develop a hydrodynamic framework to de- 
scribe energy-driven drag. The neutral modes (particle- 
hole excitations, or temperature imbalance) which me- 
diate drag in our mechanism are of a long-range char- 
acter, propagating over distances much larger than the 
inelastic mean free path I = v/j. The length scales 
relevant for our drag mechanism are the electron-lattice 
cooling length and the inhomogeneity correlation length 
£tj £dis 3> I whichever is the smallest (£ T can be as large 
as several microns, even at room temperature [T7HI5] ). 

It is instructive to compare the Feynman diagrams de- 
scribing the conventional momentum drag and our mech- 
anism (see Fig{2]) . While the characteristic momenta are 
large (~ hp) for the former contribution (Figj2^i), the lat- 
ter contribution (Fig{2|)) includes ladder diagrams repre- 
senting neutral modes propagating over distances of or- 
der £t- As a result, energy-driven drag is captured by a 
hydrodynamic framework which involves charge current 
j and heat current j q , which in the ballistic transport 
regime are related by 



j q (r) = Q(n)j, Q[n(r)] = S[n(r)]T/e, 



(2) 



where S(n) is the entropy per carrier, n(r) is the density 
profile, and e < is carrier charge. In the ballistic regime, 
using electron temperature approximation, we find 20 



Q = 



27T 2 fc|T 2 /i 

3e(p 2 + A 2 {T))' 



(3) 



where A(T) accounts for the Dirac point broadening due 
to disorder and thermal fluctuations. 

To illustrate the relation between energy and charge 
transport, we first analyze in-plane resistivity in a sin- 
gle layer. According to Eq.Q, spatial inhomogeneity 
leads to heating/cooling in the presence of uniform charge 
current (a la the Joule-Thomson process). The spa- 
tial temperature profile can be found from — VkV<5T + 



XST = — V • j q , where k is the thermal conductivity 
and XST is the electron-lattice cooling power. A tem- 
perature gradient VST drives thermopower, providing 
additional dissipation and thereby increasing resistivity. 
Onsager reciprocity combined with Eq.Q gives E(r) = 
— (Q[n(r)]/T)VffT[20]. Taking an ensemble average over 
small density fluctuations, Sp <C ksT, p, we find an in- 
crease in the in-plane resistivity, (p a p) = p^g + Ap a p, 
[a(/3) = x,y], where 



|q|<V 



(SQ(-q)SQ( q)) 
Kq 2 + X 



-q a qp- 



(4) 



Since the derivative dQ/dp peaks at p = 0, this results 
in Ap a g that peaks at CN. The temperature dependence 
estimated below is Ap oc T 2 , reminiscent of super-linear 
power laws for resistivity frequently observed at small 
doping 2 lj. A contribution of nonthermal modes to Ap 
was analyzed in Ref. |22j . 

Generalizing this analysis to two layers coupled by ver- 
tical energy transfer and accounting for correlated den- 
sity fluctuations, Eq.([T]), we find an ensemble-averaged 
drag response E 2 = P21J1, 

1 dQ dQ ^ (8p 2 (-q)Sp 1 (q L )) 



P21 



2Tk dpi d/i2 



where k = n\ + n-i is the net thermal conductivity of 
the two layers, p is the chemical potential, £ c ~ I is the 
intcrlayer cooling calculated below. 

Because the sign of the correlator in Eq.Q is positive, 
energy-driven drag has the same sign as Ap a p in Eq.Q, 
i.e. is positive at zero doping. This results in a double 
sign change along the main diagonal ri\ = ri2, as pictured 
in Fig. [l](b). The density dependence for P21 features 
a peak at zero doping (see Figj3j) which is a hallmark of 
the energy-driven drag regime. 

We note that, if the sign of the correlator (SpiSp2) 
were negative, as expected for strain-induced charge 
puddles [23] , our analysis carries through but predicts a 
negative drag at zero doping. Hence drag is a useful tool 
for probing the origin of inhomogeneity in graphene. 

We begin by studying the energy transfer between the 
electronic systems in the two layers (Figflja)). This is 
described by the Hamiltonian 



tel-el 



(6) 



where {i,j} — 1...2N index layer, and spin/valley 
degrees of freedom, Sp(r) describes the slowly vary- 
ing disorder potential, v is the Fermi velocity, 
and the electron-electron interactions are Hd-ci = 
Iiq,k,k', ! , 3 ^(q)^+ q ,ik'- q ,^k' J i/'k, ! with Vijlq) the 
interaction. 

In our analysis, we ignore the correction due to fi- 
nite layer separation d, approximating the interlayer in- 
teraction by the bare Coulomb interaction, Vy(q) = 
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'|q| 



v: 







V°e . . 

dielectric constant 



27re 2 /e|q| with e the background 
This approximation is valid when 
the lengthscale d is small compared to the screening 
length and Fermi wavelength in the layers, which is the 
case for systems of current interest [T]. The random- 
phase approximation then yields a screened interaction 
ViM) = K/l 1 - V£(ni(q,w) +n 2 (q,w))] for i,j in dif- 
ferent layers. 

We describe the energy distribution of carriers in each 
layer by a Fermi distribution at temperatures T 12 . Using 
Fermi's golden rule we can calculate the rate of energy 
exchange between the two layers (see Appendix). In the 
degenerate limit (J,i, fJ> 2 3> fc^T, we obtain the energy loss 
power for layer 1 as 
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where is the total density of states in each layer, and 
ksTo = v(2-Ke 2 / e)(vi + v%). Notably, for equal densities 
J\i does not depend on the Fermi surface size. For equal 
densities and small temperature differences between the 
layers T\ ~ T 2 , we obtain the cooling rate 

1 dj 12 _ 9C(4)fc|T 2 ln T 
C c i dT n/ih T 

where the heat capacity C c \ — tt 2 /3k 2 3 Ti'(fi) and the 
density of states v((J.) = 2fj,/ \-kTi 2 v 2 ) for the degenerate 
limit have been used. The rate 7 increases as fj, goes 
towards neutrality, but is already quite large for p: away 
from neutrality. For typical values /1 = 100 meV, T = 
300 K, the rate 7 is about 29 ps" 1 , orders of magnitude 
faster than the electron- lattice cooling rates |17H19] . 

Vertical energy transfer couples heat transport in the 
two layers, so that the layer temperatures T±, T2 obey 



-VkiV£Ti 
-Vk 2 V(5T 2 



a(£Zi 
a(ST 2 



ST 2 ) 
«i) 



AST! 
X6T 2 



-V-jg,l 

(9) 



where a = dj\ 2 / dT [see Eq. ^] and A describes electron- 
lattice cooling. We consider only a response linear in the 
applied current, j, neglecting the quadratic joule heating 
term. Inverting the coupled linear equations, we find an 
increase in temperature in layer 2, ST 2 (r), that is driven 
by current in layer 1 as 



<5T 2 (r) = - 



LiL 2 — a' 



-(ji-V)QMr),T], 



(10) 



where jq,i is the heat current, Eq.([2j, where Li — 
— Vk^V + a + X. In what follows we suppress the A 
term since electron-lattice cooling is slow. Eq.(10) then 
predicts a value for the interlayer cooling length £ c = 
y/niK 2 /[(ni + K 2 )a\. The induced temperature profile, 
<5T 2 (r), creates thermal gradients that can drive a local 
thermopower via E 2 (r) = -(Q[n 2 (r)]/T)V^T 2 . 



50'. 



021 (A) b > 



TV 







20 


AO 










10 


30 












! 


^ 


20 




1 










-10 


10 






-20 


-? 





— T=300K 

— T=100K 
— T=50K 







50 



d) 



fj, (meV) 



20 





— T=300K 




1 


Ml = V'2 


— T=200k 








— T=100K 




0.8 




— T =50 < 










ea 


0.6 














CM 


0.4 






ft. 








0.2 











-100 -50 50 

ji (meV) 




0.5 1 1.5 2 

k B T/A 



FIG. 3: (a) Total drag resistivity P21 = P21' + P21 
chemical potentials in the two layers, evaluated from Eq. allh 
and Eq.|5| at T = 100 K, producing a peak at /Lti,2 = (see 
text for parameter values used). (b,c) Slices fii — [1,2 and fii = 
—fj,2 at different temperatures. Note a three-peak structure 
in slice (b) and two sign changes close to CN in (c). (d) 
Temperature dependence of the peak at /Lti,2 = 0. 

Spatial fluctuations in thermopower are governed by 
density fluctuations via Eq.(lO). In particular, close 
to neutrality the local thermopower will exhibit regions 
of both positive and negative sign, leading to a spa- 
tial pattern of the drag resistivity. As discussed above, 
the correlations between SfXi and Sfx 2 , Eq.Q, lead to a 
nonzero ensemble-averaged drag resistivity. In the limit 
SfJ,i,2 < k B T,m, 2 we write Q l (r) = (Q l (r)) + ^5^{r). 
Passing to Fourier harmonics via (Sfj,i(r)5u 2 (r')) = 
E q e <q(r_r ' ) <^i(-q)^2(q)>, we obtain Eq.jHj). 

The fact that fluctuating local thermopower, exhibit- 
ing both positive and negative signs, does not average to 
zero is surprising. This happens because the inhomogene- 
ity in heat current and thermopower arise from the same 
source: electron-hole puddles. Energy-driven drag resem- 
bles mutual drag described by Laikhtman and Solomon 
[2"4~] in semiconducting heterostructures where doping at 
contacts produced a similar correlation between Peltier 
heating/cooling and thermopower. Energy-driven drag 
in Graphene differs from Ref.[23] in that density inho- 
mogeneity is intrinsic, occurs throughout the sample (not 
just at the contacts), and on a far smaller scale. 

To see how energy-driven drag, Eq. ([5| , affects the total 
experimentally measured drag we need to account for 
momentum drag, . We use a model that captures 
the main qualitative features of momentum drag: 
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This expression, with the prefactor p 21 = 
— 1.4a 2 /(27r?7 2 ), the effective interaction strength 
a = 0.05, and 77 = 6.25, was obtained by fitting the 
functional dependence derived in Ref. [7] in the doping 
region -10 < p/k B T < 10. 

Combining this with p^i m Eq.dHb, we obtain the to- 
tal drag p 2 °^ = + /4i plotted in Figji] Here we 
have used k = (p 2 + A 2 (T))/hT [25], Eq.Q and assumed 
Gaussian correlations with average square density fluctu- 
ations {8 ii 2 ) fa 25meV 2 and £ dis = 100 nm [HI [IB]. We 
note that the functional form of the correlator in Eq.Q 
does not impact the qualitative behavior. The obtained 
values of total drag are compatible with measured drag 
resistivities reported in Refs.[5l HP]. 

The density dependence of total drag plotted in Fig[3] 
(a) can be used to distinguish the two drag mechanisms 
in experiments. Namely, the peak at zero doping is due 
to energy-driven drag. On the slice jj>x = —p>2 (black 
dashed line) this peak is surrounded by two peaks domi- 
nated by the momentum contribution [Fig[3jb)]. On the 
slice /ii = fj,2 (yellow dashed line) the two mechanisms 
produce contributions of opposite sign, resulting in a dou- 
ble sign change [Fig{3|c)]. This provides a clear means 
of discerning the energy-driven regime. 

The temperature dependence can be estimated as fol- 
lows. At very low T such that £,£ c > £dis 7 the sum 
in Eq.([5j is cut at l/£, giving P21 oc T 8 . At not too 
low T such that £ c ,£ < £dis, the sum in Eq.([5]) yields 
£ q (<W-q)<Wq)) = (Spi(r)Sp 2 {r')) r =r'- This gives a 
non-monotonic T dependence [see Fig. [3] (d)] 

P21 oc — — 3^/i 1 (r)^ 2 (r'))r=r'- (12) 

(A 2 + ry(fc B T) 2 ) 3 

A similar non-monotonic T dependence arises for in- 
plane resistivity A.p a p. Interestingly, a peak in drag re- 
sistance at CN with non-monotonic temperature depen- 
dence was recently reported in Ref.|10j. 

The above analysis can be easily extended to describe 
the diffusive limit where the elastic mean free path is 
shorter than the inelastic mean free path, £' < I. Our hy- 
drodynamic approach remains valid in this regime, with 
the quantity Q = sT where s is the Seebeck coefficient. 
The energy-driven drag is still given by Eq.([5]), with s 
and k described by the Mott and Wiedemann-Franz re- 
lations: 

s -^k 2 T— e 2 K -—k 2 Ta (IS) 
S ~ Se kBl dp ' 6 K ~ 3 ksia ' {U} 

where a is the electrical conductivity. Taking a to vary 
linearly with carrier density, we find Q in the disordered 
limit takes on the same qualitative form as Eq. ^ in the 
clean limit. As a result, the qualitative features of energy- 
driven drag in the clean limit also appear in the disor- 
dered limit: namely, the avoided crossing of nodal lines 
(Fig. [31(a)), a peak at zero doping, double sign reversal 



along the diagonal ni = n 2 and a three-peak structure 
along the diagonal n\ — —n 2 (Fig. [3](b,c)). The T de- 
pendence of P21 is qualitatively similar in the diffusive 
and ballistic regimes. However, since the Wiedemann- 
Franz relation gives n oc T (as opposed to k oc 1/T in 
the ballistic regime), we find P21 oc T 2 at lowest T and 
P21 oc T~ 4 at higher T > p,, A, as shown in Fig. [3] (d). 

In summary, vertical energy transfer in graphene het- 
erostructures has strong impact on lateral charge trans- 
port in the Coulomb drag regime, dominating the drag 
response at low doping. Drag measurements thus afford 
a unique probe of energy transfer at the nanoscale, a fun- 
damental process which is not easily amenable to more 
conventional techniques such as calorimetry, and is key 
for the physics of strong interactions that occur near neu- 
trality. 
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APPENDIX A: HEAT CURRENT AND 
ONSAGER RECIPROCITY 



Here we use Onsager reciprocity to obtain ther- 
mopower induced by a temperature imbalance in a Fermi 
gas. This can be done in a general form applicable in 
both the ballistic regime, when the mean free path is 
dominated by electron-electron scattering, and in the dif- 
fusive regime, when the mean free path is dominated by 
elastic scattering by disorder. 

We start with recalling Onsager reciprocity for the heat 
and charge transport. Given charge current j and heat 
current j 9 described by 



-j = £ll^V/Z + £l2V^ 
jq = L 2 l^VfJ, + L 22 \l} 



T 



T 



(14) 
(15) 



the cross-couplings obey the Onsager relation L\ 2 = Lai- 
Next, we consider the heat current, which can be ob- 
tained from the heat transport equation 

d t (C cl T) - VkVT = -j • VF[n(r),T], (16) 

where C c \ is the electron heat capacity, F[n(r),T) is a 
function of the carrier density n(r) to be determined 
later, T is the temperature and k is the thermal con- 
ductivity. From Eq.(16) we find that the heat current 
is 

j q = -KVT+jF[n(r),T] = -nVT - aF[n(r), T]V/i, 

(17) 

which gives L 2 \ — — TaF[n(r), T] in Eq.([l4|. Here we 
have used V • j = and j = — erV/i, where a = L\{T. 
Using the Onsager relation, the thermopower e.m.f. in- 
duced by the temperature gradient equals 



E 



F[n(r),T] 



VT. 



(18) 



This result has general validity irrespective of the trans- 
port mechanism specifics, which are manifested through 
the form of F[n(r),T]. This quantity equals the See- 



beck coefficient s given by Eq.( 13 1 in the diffusive regime, 
whereas in the ballistic regime it is given by Eq.(|3]). 

The functional form of F[n(r),T] can be obtained by 
considering the kinetic equation (at steady state) 



eE- V p n(p,r) = h + 1 2 



(19) 



where e < is the carrier charge. Here we write the colli- 
sion integral as a sum of momentum non-conserving and 
momentum conserving parts, I\ + I 2 , corresponding to 
disorder scattering and electron-electron scattering, re- 
spectively. Heat and change current can be expressed 
through a steady-states deviation from the equilibrium 
Fermi distribution, Sn — n — no, as follows 

jq = X^ 6 * ~ A t )v P ,^n i (p,r), j = v P!l <Sn 4 (p, r). 



p. 1 



(20) 

Here i labels the conduction and valence band states and 
v p,i = dei/dp and ej = ±u|p|. Below we consider the dif- 
fusive and ballistic regimes. In the first case the mean free 
path is dominated by clastic momentum non-conserving 
scattering (/1), in the second regime the mean free path is 
dominated by inelastic momentum-conserving scattering 
(h). 

In the diffusive regime, neglecting I 2 and using the 
relaxation time approximation for I\, we find 



= eE 



v p T(e)dn F / de, 



(21) 



where r describes elastic scattering by impurities. This 
gives the standard expressions for Seebeck and Pelteir 



coefficient described in Eq.(13l of the main text, so that 
F = sT. 

In the ballistic regime, the fastest scattering mecha- 
nism comes from the (total) momentum conserving pro- 
cess of electron-electron scattering. As a result, we will 
neglect all other terms apart from I 2 in Eq. [19] and look 
for distributions, n(p,r), that give a non-zero particle 
flow. At a nonzero total current, the non-equilibrium 
distribution can be written as np(e — P ■ u), where the 
term p • u describes the change due to particle flow. This 
allows us to write Sn as 

drip 



up 



de 



(22) 



where the Fermi distribution, np , has a temperature that 
may depend on the flow and position. Using Eq. [22] and 
summing over both conduction and valence bands (where 
e c = v\p\ and e v — — v\p\, c and v refer to conduction and 



valence bands respectively) in Eq. 20 we obtain 

. drip 



u 

~2~ 
— ue 

~2~ 



dev(e)e{e — ji)- 



de 



dev(e)e 



dnp 



(23) 
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where v(e) — 2\e\ / {nfi 2 v 2 ) is the total density of states 
and we take into account that v p ■ p = e with e posi- 
tive (negative) for the conduction (valence) band. Since 
the function ev(e) oc e|e|, the integral cannot be evalu- 
ated for arbitrary ratio ksT/p. Instead, we analyze the 
degenerate case, p. fc^T, in which case we use the 
Sommerfeld expansion to obtain 

j=| M KM)u, j q = 2C(2)^)(fc B T) 2 u (24) 

where we have used the identity / °° = C(2) = 

7r 2 /6. Comparing both expressions to eliminate u we 
obtain 



e/i 



(25) 



The singularity at fi — is smeared by broadening of the 
Dirac point due to disorder and thermal fluctuations. We 
can account for smearing via 



. 4C(2) (k B T) 2 p 
Jq e /i 2 + A 2 (T) 



j, A 2 (T) = A 2 + ^(fc B r) 2 (26) 



where A is the Dirac point width parameter. This yields 
F[n(r), T] = Q, giving Eq.Q of the main text. 



APPENDIX B: INTERLAYER ENERGY 
RELAXATION 

In this Appendix, we derive the interlayer cooling 
power for the heterostructure system defined by the 
Hamiltonian in Eq.([6]). Using Fermi's golden rule, the 
scattering rate is 



k x ',ki 



2ttN x ^ 

k 2 ,k 2 ' 

q 



F kaMa Wf(e lsa )[l-f(e k , l )]6 e 8 1 6 2 



(27) 

where S e = S(e kl > Ml + £k 2 ',k 2 ), and 5 1 = S kl > Ml+cl and 
$2 = ^k 2 '.k 2 -q- Here N = 4 is the number of spin/valley 
flavors, {1, 2} denote the different layers, ek',k = £k' — £k, 
and Fk.k' = |(k'a|k/3)| 2 is the coherence factor (a, f3 
label states in the electron and hole Dirac cones). V q 
is the screened inter-layer Coulomb interaction described 
below. The energy-loss power is 



J = -N 



E W' klSkl (e ki -e kl )/(c kl )[l-/(c ki )]F klikl , 



ki,ki 

(28) 

We can simplify the evaluation of these sums by writing 
S{e kl >, kl +e Wik2 ) = SZ d^6(e kl > Ml - + e k2 ',k 2 ) 

and using the identity / s (e k )[l - / s (e k ')] = (Js(ek) - 
/ s (e k ')) x {N s {e k , tk ) + 1) where N(u) = l/(e^ T - 1) 
is the Bose function taken at the electron temperature (of 
that particular layer) and s = {1,2} denotes the layers. 



Using the quantities 

X"(q,w) = N V J Fk,k+q(/ s (ek)-/ s (ek+q))'5(ek+q-ek-w). 



(29) 

where x"(q, w) = ilmll(q, us) is the imagi- 
nary part of the susceptibility and fl(q, lo) = 

N Ek *k,k +q {^l-l k -?o is the Polarization op- 
erator. Using these we can re-write the energy loss 
power as 

J=lJ dww(JV a (u)-JVi(w))2|^ q |M(q.w)X2(q,w) 

(30) 

where we have noted that N(—uj) = —(1 + N(u>)) and 
used the fact that x"(q 5 w ) = — x"( — Qi ~ w ) so as to only 
keep the odd part of the product N 2 (— uj)(Ni(uj) + 1). 

The energy transfer between the two layers is depen- 
dent on the Coulomb interaction Vq between the layers. 
We treat V q by accounting for polarization in both layers 
and screening in the RPA approximation. The RPA- 
screened coulomb interaction is 



Fq= l-V?(n 1 (q,i+n 2 (q,a;)) ! 



VZ = — (3D 
H eq 



where s is the background dielectric constant. Here we ig- 
nored the correction due to finite interlayer spacing <i, ap- 
proximating the interlayer interaction U q = -^- e - d \<i\ ~ 
Uq . This approximation is valid when the layer separa- 
tion d is small compared to the screening length in the 
layers, dq < 1, q = -Ili(q = 0) - n 2 (q = 0). We will 
be interested in the regime when d is small compared to 
the Fermi wavelength in the layers, for which the above 
approximation is adequate. 

In the degenerate limit, /i ^ fcsT, the polarization is 
given by 



U(q,u) = 1 



yj{uj + i0) 2 - q 2 v 2 



(32) 



In the limit of q <C qo where qo = {2-iTe 2 j e){y\+v 2 \ where 
Vi 2 is the density of states at the Fermi level in each 
layer, we can write \ V q \ = l/|II 1 (q, w)+n 2 (q, w)\. Noting 
that in the degenerate limit, intra-band transitions are 
the dominant processes we approximate x"(q, u>) by the 
imaginary part of Eq.(32). This allows us to write the 
cooling power as 



J12 = \ J duiuj(N 2 (Ls) - Ntiu)) J2 



V\V 2 UJ 



|q|>w/u 



(i>i + v 2 ) 2 q 2 v 2 



(33) 

where To = vqo/ks- Using the idenity J_ uj 3 du>N(uj) — 
2 ■ 3!C(4)(fcsT) 4 , and recalling that the ultra-violet cutoff 
arises from screening via q -C go we obtain Eq.([7|. 



